NASA Technical Memorandum 88803 


A Mathematical Approach for Evaluating 
Nickel-Hydrogen Cells 


(NASA-TH-88803) A M ATBEH ATICA1 
evaluating nickel-hydrogen cells 
18 P 


APPROACH FOR 
(NASA) 

CSC L 10C 


G3/25 


H 86-3 16 80 


Onclas 

43502 


Harold F. Leibecki 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

1986 Fall Annual Meeting of the AIChE 
Miami, Florida, November 2-7, 1986 





E- 3114 
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ABSTRACT 

A mathematical equation Is presented which gives a quantitative 
relationship between time-voltage discharge curves, when a cell's ampere-hour 
capacity Is determined at a constant discharge current. In particular the 
equation quantifies the Initial exponential voltage decay; the rate of voltage 
decay; the overall voltage shift of the curve and the total capacity of the 
cell at the given discharge current. The results of 12 nickel-hydrogen boiler 
plate cells cycled to 80 percent depth-of-dlscharge (DOD) are discussed In 
association with these equations. 

INTRODUCTION 

The shape of a cell's discharge determination (100 percent DOD) Is the 
result of complex Interactions of several rate limiting processes. These 
processes may be physical or chemical In nature and changeable during the life 
of the cell. The anode, cathode or electrolyte phase may be responsible 
Individually or In combination for the shape of the cell's discharge curve. 
Although specific tests have been designed to lllucldate certain Isolated 
characteristic of the discharge curve, the customary practice of displaying 
complete discharge curves Is still used to show changes In cell behavior 
during life and differences between cells. 

The changes associated with discharge curves give Insight to the 
electrochemist as to the gross properties of the cell. However; when many 
cells are being Investigated the data processing and trend analysis procedure 
become unmanageable without the use of mathematical techniques. By the use of 
an equation that describes the various characteristics of, the discharge shape 



It Is possible to assign definitive attributes to the discharge curve. This 
then reduces the amount of Information to be assimilated to a manageable 
amount. This Information now being In numerical form allows the use of other 
mathematical procedures to further lllucldate the changes In the curves. 

Two related mathematical equations are used to compare changes occurring 
during cycling as well as differences between cells for 12 nickel-hydrogen 
cells. The numerical values from the equations are further analyzed with 
respect to cycle life. 

EXPERIMENTAL DATA 

Nineteen nickel-hydrogen boiler plate cells were constructed and life 
cycled by Hughes Research Laboratories (Ref. 1). These cells were constructed 
from nickel plaque having three different bend strengths and three different 
pore sizes. Nickel electrodes were prepared from these plaques having three 
different loadings of active material. Table I shows the cell number and Its 
respective design characteristics. Initial performance tests Included 12 
consecutive capacity measurements to characterize the cells. To facilitate 
experimental cycling the cells were divided Into three rated capacity groups 
of 2.7, 3.0, and 3.3 A-hr. This grouping allowed for cycling all cells using 
a constant time for charge and discharge with three different currents, 
dependent upon the rated capacity. 

The cells were life cycled at 23 °C to 80 percent depth-of-dlscharge (DOD) 
of the rated capacities using a 45 min cycle regime. The cycling regime 
consisted of a 2.74 C rate discharge for 17.5 min and a 1.92 C rate charge for 
27.5 min. Thus, as near as possible normalizing the test conditions. The life 
cycling regime was Interrupted periodically for capacity measurements after 
approximately every 1600 cycles. Capacity measurement data at cycles 1700, 
3100, 5000, and 6100 were available as strip chart recordings of voltage-time 
curves for the 1.37 C discharge rate. The data was transformed Into numerical 
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voltage-time points by digitizing the strip chart using a Hewlett-Packard 9111A 
'-.Graphics Tablet. The accuracy of this process was found to be ±2 mV and 
±0.1 min. This results In a maximum uncertainty of ±0.002 A-hr. 

The conversion from voltage-time to voltage-OOD was performed so as to 

}:• • 

further normalize the variation In cell capacity. The pre-cycling capacity at 
1.37 C rate discharge was used In this normalization process since the test 
data was also at the 1.37 C rate. Equation (1) describe the conversion. 

n (time (min) * current (All 

000 = 60 mln/hr * cap @1.37 C (A-hr) (1) 

EQUATION 

It was found (Ref. 2) that Eq. (2) describes the typical s-shaped 
discharge profile 

E • K - shr ' Ln [^(ij * A o * Exp( - A i * 000 » < 2 > 

where 

X D0D/A 

A Capacity of the cell at the specified discharge conditions (D0D) 

Ln(C) Is proportional to the slope of the discharge curve (volts ^ 

Ag Maximum height of the voltage drop at the start of discharge (volts). 

A.j the voltage drop duration at the start of discharge (00D -1 ) 

K Intercept of the voltage plateau with the voltage axis taken at 

D0D = Exp(-1 ) (volts) 

These relationships are shown In Fig. 1. 

Equation (2) can be expanded to describe the discharge profile of a two 
plateau curve as shown In Eq. (3). 
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where: All parameters have the same significance as In Eq. (1). The (') 

parameters are associated with the second plateau. 

The digitized voltage-DOD data was fit to Eq. (3) by using a nonlinear 
regression procedure on a Hewlett-Packard 9845C desk top computer. The 
regression program follows the Marquardt procedure to obtain estimates for the 
parameter. 

If only one plateau Is present In the voltage-DOD curve Eq. (2) can be 
used. However, all cells exhibited two plateau discharge characteristics, to 
various degrees between open circuit and 0.4 V. Therefore, Eq. (3) was used. 
Figure 2(a) shows the result of boiler plate cell number 12 at the 5000 cycle 
capacity check which has a very small second plateau (A'-A = 0.068 D0D). 

Figure 2(b) shows the result of boiler plate cell number 14 also at the 5000 
cycle capacity check which has a moderate well formed second plateau 
(A'-A = 0.152 D0D). Figure 2(c) Is a plot of the difference between the 
experimental potential and the results of Eq. (3) for each D0D point digitized 
for cell number 14 at the 5000 cycle capacity check. As can be seen the 
equation provides a good representation of the data. The average of the 
potential differences for 23 data points Is -2.5 mV with a standard deviation 
of 6.5 mV. 

ANALYSIS OF DATA 

The cycle life data as shown In Table I varies from 2330 to 12 960 cycles. 
This decade of variance, however, cannot be attributed to any of the 
experimental design characteristics (bend strength, pore size or loading) since 
no statistically significant (1-a = 0.9) difference could be found using 
analysis of variance techniques. With a cycle life variation of over 10 000 
cycles and no effect due to experimental design It Is unlikely that cell cycle 
life Is 6185 with a standard deviation ±3054 cycles. It Is therefore necessary 
to look elsewhere for the cause of variation. 
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Data (E-DOD) was available for 12 of the 19 cells at capacity 
determination cycles 1700, 3100, 5000, and 6100. 1 Boiler plate cells 
number 23, 24, and 25 were early failures (approximately 2300 cycles) and 
therefore were not used. Boiler plate cells number 16, 18, 19, and 26 were 
available for only one or two capacity determination cycles and were also 
excluded. All 12 boiler plate cells had capacity check data for capacity 
measurements at cycles 1700, 3100, and 5000; however, six cells had cycle lives 
long enough to have a capacity measurement at cycle 6100. All 42 capacity 
determinations were fitted to Eq. (3) and the parameters of the equation used 
as a means of evaluating cycle life. 

l 

H.S. Urn (Ref. 1) attempted to fit the cycle life of cells to the percent 
D0D (based on the theoretical capacity). Although a trend was observed the 
amount of scatter was too large to draw conclusive results. The use of the 
theoretical D0D would have to assume that cells would behave Identically during 
their early lifetime l.e., gain or lose capacity at the same rates during 
formation. Also, reported by Llm was the appearance of the second plateau that 
varied In magnitude In a unrecognizable manner during life. 

From the Information obtained by fitting the data to Eq. (3) It Is 
apparent that A' varies during life. This Indicates that the capacity of the 
cells are neither constant nor decreasing unlnformally during their cycle 
lives. For cells with short cycle lives, A' tends to continually decrease 
over life. For cells of Intermediate life (4000 to 5000 cycles) the trend Is 
towards no change In the first two capacity checks (1700 & 3100) and then a 
decrease In capacity at the 5000 cycle/check. For cells with long life, (6000 
and more cycles) the trend Is for an Increase In A' and then a decrease. 


^No capacity determinations were performed beyond the cells cycle life. 
If a cell failed at cycle 3000 capacity determination at 3100 was the last 
measurement. 
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In addition to the TO equation parameters obtained from the nonlinear 
data fit to Eq. (3) two additional quantities were used. These additional 
quantities are the DOO's to which the cells were actually cycled and the ratio 
of A-hr at the capacity determination cycle to the cell's theoretical A-hr 
capacity (utilization). Both of these values are based upon the parameter A' 
from Eq. (3) and are shown In Eqs. (4) and (5). 


Ampere-hour ratio (Ut) . ^ ~ etlcIl^lacUy 1 
cycling D00 (Cy 000) - f r 'i ^firfl3Tc) 


(4) 

(5) 


where capacity (1.37 C) = pre cycling capacity at 1.37 C discharge rate. 

Figure 3 shows how the capacity of the cells have Increased above the 
theoretical capacity at capacity determination cycle 1700 and then decreases 
as cycle life progresses. This Increase Is not phenomenal or unexpected. It 
Is a well observed fact that nickel electrodes Increase capacity during Initial 
cycling (Refs. 3 and 4). This change In capacity alters the actual 000 of a 
cell during cycling and Is shown In Fig. 4. As can be seen from Fig. 4 the 
Increase In capacity has lowered the actual D0D of cycling to about 65 percent 
at capacity determination cycle 1700. From the capacity determinations at 
cycles 3100 and 5000, the variation In cell cycling 000 Increases to three 
times that at cycle 1700 (0.1 to 0.311). Since cycle life Is strongly 
dependent upon 000 this variation In actual 000 phenomena should effect cycle 
life. 


In order to determine the Impact of the changing discharge curve on cycle 
life the 10 parameters from Eq. (3) along with the two quantities describing 
the A-hr ratio (Ut) and cycling 000 (Cy 00D) were regressed to cycle life using 
a stepwise forward linear regression procedure (Ref. 5). Equations (6) to (9) 
show the results of the linear regression for each of the capacity 
determination cycles. 
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Capacity determination at cycle 1700: 

Life = 102 902 - 16 840 * A - 36 903 * Ut - 63 969 * Cy D00 (6) 

Capacity determination at cycle 3100: 

Life = 123 073 - 33 931 * A - 34 298 * Ut - 76 707 * Cy DOD (7) 

Capacity determination at cycle 5000: 

Life = 61 155 - 31 927 * A Q - 23 699 * Ut - 37 461 * Cy 000 (8) 

Capacity determination at cycle 6100: 

Life = 29 558 - 54 821 * A fl - 17 666 * Cy DOD (9) 


As can be seen from Eqs. (6) to (9) actual DOD of cycling Is the major 
contributor In all but the capacity determination at cycle 6100. Figures 5(a) 
to (d) show plots of how well Eqs. (6) to (9) predict cycle life. Although 
Eqs. (6) to (9) do not produce a perfect fit to cycle life, as portrayed by 
the doted line, 88, 90, 90, and 99 percent of the variability In cycle life 
can be accounted for by Eqs. (6) to (9) respectively. Table II lists the 
difference between the actual and predicted cycle lives for the four capacity 
determination cycles. With the exception of boiler plate cell number 20 the 
difference Is acceptable as a predictor of cycle life. Table III lists the 
percentage difference between the actual and predicted cycle lives. Boiler 
plate cell number 20 Is shown to be In error by over 50 percent for all three 
capacity determination cycles. The mean of the percent residuals with boiler 
plate cell number 20 excluded Is only -2, -3, -1, and -2 percent for capacity 
determinations at cycles 1700, 3100, 5000, and 6100 respectively. The reason 
boiler plate cell number 20 Is so much different than the other Is not known. 

Examination of Eqs. (6) to (9) show some rather Interesting trends. 
Parameter A which Is a measure of the capacity available at the higher 
potential Is a predictor early In life of cycle life while parameter Ag 
which Is a measure of the potential drop from open-circuit Is not. This trend 
reverses Itself after capacity determination cycle 3100. The ampere-hour ratio 
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(utilization) decreases In Importance as cycle life Increases droplng to zero 
at capacity determination cycle 6100. The actual cycling DOD Is the only 
contributor for predicting cycle life needed In all capacity determination 
cycles. 

CONCLUSIONS 

Although the response variable, cycle life, showed no significant 
dependence upon the experimental variables bend strength, pore size or loading 
level, the discharge curves of capacity measurements during life, are related 
to cycle life. 

It was shown that: 

1. Discharge curves can be represented by an equation even when there are 
two plateaus. 

2. Parameters from the equation describing the discharge curve are related 
to cycle life. 

3. Cell cycling depth-of-dlscharge and Ah ratio (utilization) are very 
Important characteristics which can be obtained from the equation describing 
the discharge curve. 

4. The capacity obtained during cell characterization Is not a reliable 
Indicator of cell capacity during life testing. 

5. The capacity of the higher potential plateau Is Important In predicting 
cycle life early In cycle testing. 

6. The exponential potential drop Is Important In predicting cycle life In 
the latter part of cycle testing. 

The single most Important type of information needed to understand this 
set of cells, and probably all nickel-hydrogen cells, are the processes that 
lead to the Increase of cell capacity over an extended period of cycling. 
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TABLE I. - NICKEL ELECTRODE CHARACTERIZATION PLAN 
AND CYCLE LIFE 


Loading, 

g/cc 

Bend strength, ps 1 

400 

550 

700 

Pore size, y 

10 

12 

14 

10 

12 

10 

12 

1 .40 


m 

25 

2330 

22 

3349 

17 

3664 


m 

1.55 

69 

13 
1 2960 

24 

2580 

B 

m 

o 

9 

8936 

1 .70 

a 14 

b 4790 

11 

5167 

23 

2340 

20 

4259 

m 


m 


a F1rst line denotes cell number. 
^Second line denotes cycle life. 


TABLE II. - DIFFERENCE BETWEEN ACTUAL AND 
PREDICTED CYCLE LIFE 


Cycle Difference, Actual-Predicted 


Cell 

number 

Capacity determination cycle 

Actual 

cycle 

life 

1700 

3100 

5000 

6100 

8 

-763 

-169 

923 

-29 

9 572 

9 

43 

-94 

-678 

-64 

8 936 


-992 

-350 

783 

— 

3 729 


746 

20 

-669 

-182 

5 167 

B B 

468 

-175 

-130 

— 

3 450 

■ B 

1531 

337 

820 

-57 

12 960 

19 

-30 

213 

-477 

— 

4 790 

■ 1- 

-103 

1345 

134 

-419 

7 919 

mm 

255 

318 

-417 

— 

3 664 

■fl 

-2377 

-2414 

-2172 

— 

4 259 

KB 

1300 

1341 

1067 

-86 

6 426 

22 

197 

-372 

818 

— 

3 349 




















































TABLE III. - PERCENTAGE ERROR OF PREDICTED 
CYCLE LIFE 

Percent Error TOO*, Actual-Predicted/Actual 


Cell 

number 

Capacity determination cycle 

Actual 

cycle 

life 

1700 

7100 

5000 

6100 

8 

-8 


9.6 

-0.3 

9 572 

9 



-7.6 

-0.7 > 

8 936 

10 



21 

— 

3 729 

11 



-12.9 

-3.5 

5 167 

12 

13.6 

-5.1 

-3.8 

— 

3 450 

13 

11.8 


6.3 


12 960 

14 

-0.6 


-10 

— 

4 790 

15 

-1.3 


1.7 

-5.3 

7 919 

17 

7 


-11.4 

— 

3 664 

20 

-55.8 


-51 

— 

4 259 

21 

20.2 


16.6 

-1.3 

6 426 

22 

5.9 

-11.1 

24.4 

— 

3 349 

Mean 3 

-2 

-3 

-1 

-2 



a Cell number 20 not used In mean calculation. 
























VOLTAGE DIFFERENCE (mv) 








1.2 



1700 3100 5000 6100 

CAPACITY DETERMINATION CYCLE 


FIG. 3 Ratio of cell actual capacity to theoretical capacity. 
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FIG. 4 Cell actual cycling DOD. Design DOD = 0.8 
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FIG. 5a Predicted cycle life from 1700 cycle capacity 
determination data. 
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FIG. 5b Predicted cycle life from 3100 cycle capacity 
determination data. 




PREDICTED CYCLE LIFE (10 ~ 3 ) PREDICTED CYCLE LIFE (10~ 3 ) 



FIG. 5c Predicted cycle life from 5000 cycle capacity 
determination data. 
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Predicted cycle life from 6100 cycle capacity 
determination data. 
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FIG. 5d 





2. Government Accession No. 


1. Report No. 

NASA TM-88803 


4. Title and Subtitle 


A Mathematical Approach for Evaluating 
Nickel -Hydrogen Cells 


7. Authors) 

Harold F. Lelbeckl 


3. Recipient's Catalog No. 


5. Report Date 


6. Performing Organization Code 

506-41-21 


8. Performing Organization Report No. 

E— 31 1 4 


to. Work Unit No. 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


1 1 . Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Memorandum 


14. Sponsoring Agency Code 


IS. Supplementary Notes 

Prepared for the 1986 Fall Annual Meeting of the AIChE, Miami, Florida, November 
2-7, 1986. 


16. Abstract 

A mathematical equation Is presented which gives a quantitative relationship 
between time-voltage discharge curves, when a cell's ampere-hour capacity Is 
determined at a constant discharge current. In particular the equation quanti- 
fies the Initial exponential voltage decay; the rate of voltage decay; the over- 
all voltage shift of the curve and the total capacity of the cell at the given 
discharge current. The results of 12 nickel-hydrogen boiler plate cells cycled 
to 80 percent depth-of-dlscharge (OOD) are discussed In association with these 
equations. 


17. Key Words (Suggested by Authors)) 

Batteries 
Fuel cells 


18. Distribution Statement 


Unclassified - unlimited 
STAR Category 25 


19. Security Classlf. (of this report 

Unclasslf le 



20. Security Classlf. (of this page) 

Unclassified 


21. No. of pages 122. Price* 


*For sale by the National Technical Inlormaiion Service. Springfield, Virginia 22161 


















National Aeronautics and second class mail 

Space Administration 

Lewis Research Center address correction requested 

Cleveland. Ohio 44135 

Official Business 
Penalty lor Private Use $300 



Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 4 


I 




I 



